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Abstract 

We preliminarily investigate the impact of the Kuiper Belt Objects (KBOs) and of the asteroid ring on some proposed high- 
precision tests of Newtonian and post-Newtonian gravity to be performed in the Solar System by means of spacecraft in heliocentric 
^ ~ 1 AU orbits and accurate orbit determination of some of the inner planets. It turns out that the Classical KBOSs (CKBOS), 
^— ^ which amount to ~ 70% of the observed population of Trans-Neptunian bodies, induce a systematic secular error of about 1 m 
after one year in the transverse direction T of the orbit of a test particle orbiting at 1 AU from the Sun. For Mercury the ratios 
f— «) of the secular perihelion precessions induced by CKBOs to the ones induced by the general relativity and the solar oblateness J2 
f — , amount to 6 x 10~ 7 and 8 x 10~ 4 , respectively. The secular transverse perturbation induced on a ~ 1 AU orbit by the asteroid 
f*"*) ring, which globally accounts for the action of the minor asteroids whose mass is about 5 x 10~ 10 solar masses, is 10 m yr -1 ; the 
""^ bias on the relativistic and J2 Mercury perihelion precessions is 6.1 x 10~ 6 and 1 x 10 -2 , respectively. Given the very ambitious 
^-igoals of many expensive and complex missions aimed to testing gravitational theories to unprecedented levels of accuracy, these 
I notes may suggest further and more accurate investigations of such sources of potentially insidious systematic bias. 
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1. Introduction 

The goal of this letter is to draw attention on the poten- 
tially corrupting effect of the Kuiper Belt Objects (KBOs) 
and of the ring of the minor asteroids on many proposed 
high-precision tests of Newtonian and post-Newtonian 
gravity to be performed in a more or less near future with 
spacecraft moving along heliocentric orbits in the inner 
regions of the Solar System. 

2. The PPN parameters (3 and 7 and their current 
experimental bounds 

In the Parameterized Post-Newtonian (PPN) formalism 
(Will, 1993), developed by Nordtvedt and Will to make 
easier the comparison of metric theories of gravity with 
each other and with experiment, the Eddington-Robcrtson- 
Schiff parameters f3 and 7 describe how much non-linearity 
is present in the superposition law for gravity and how much 
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spatial curvature is produced by a unit mass, respectively; 
in general relativity (3 = 7 = 1 . For a recent overview of the 
measurements of such PPN parameters see (Will, 2006). 

The most accurate and clean determinations of 7 come 
from the effects involving the propagation of electromag- 
netic signals: in their PPN expressions only 7 is present, 
at first order. For the Shapiro time delay (Shapiro, 1964), 
recently tested with the Doppler tracking of the Cassini 
spacecraft, the most accurate result is 7 — 1 = (2.1 ± 2.3) x 
10~ 5 (Bcrtotti ct al., 2003). With (3 the situation is less fa- 
vorable because it enters the post-Newtonian equations of 
motion of test particles in conjunction with 7. For exam- 
ple, in the post-Newtonian expression of the well known 
Einstein precession of the pericentre u> of a test particle 
(Einstein, 1915) the combination v = (2 + 27 — /3)/3 is 
present. E.V. Pitjeva recently processed more than 317 000 
Solar System observations (1913-2003) for the construction 
of the Ephemerides of Planets and the Moon EPM2004 de- 
termining, among other things, f3 and 7 with an accuracy of 
10 -4 (Pitjeva, 2005). The same limit for j3 can also be ob- 
tained by combining the Cassini result for 7 with the Lunar 



Laser Ranging test of the Nordtvedt effect 1 (Nordtvedt, 
1968a, b), expressed in terms of r\ — 4/3 — 7 — 3: indeed, 
from 5i] = 4.5 x 10~ 4 (Bertolami et al., 2006), it follows 
5(3 = 1.3 x 10~ 4 . 

3. Proposed space-based missions for 
unprecedented tests of post-Newtonian gravity 

Various theoretical models involving scalar-tensor sce- 
narios (Damour and Nordtvedt, 1993a, b; Damour and 
Esposito-Farese, 1996a, b; Damour ct al., 2002a, b) predicts 
deviations from unity for (3 and 7 at a 10 -6 — 10~ 7 level, 
so that it is very important to push the accuracy of such 
PPN tests towards this demanding accuracy. As a conse- 
quence, many space-based dedicated missions have been 
so far proposed; they share the same operational scenario, 
i.e. the inner regions of the Solar System from 0.38 AU to 
about 1 AU. Below we briefly review some of them. 

The LATOR (Laser Astrometric Test Of Relativity) mis- 
sion (Turyshev et al., 2006) has, among its goals, the de- 
termination of 7 with an accuracy of 10 -9 from the first 
order and a direct and independent measurement of (3 at a 
~ 0.01% level via the second-order gravity-induced deflec- 
tion of light. LATOR is based on the use of the Interna- 
tional Space Station (ISS) and of two small spacecraft to 
be placed in heliocentric ss 1 AU 3:2 orbits resonant with 
that of the Earth. The nominal mission lifetime amounts 
to 22 months. 

The most accurate method available for determin- 
ing planetary orbits is based on accurate range mea- 
surements from the Earth to a spacecraft orbiting an- 
other planet or to a lander on the planet's surface. The 
ESA BepiColombo mission (http://sci.esa.int/science- 
e/www/area/index.cfm?fareaid=30), which should be 
launched in 2012, places itself in such a framework aiming 
at measuring (3 and 7 at a 10 -6 level (Milani et al., 2002) 
from an accurate determination of the orbits of Mercury 
and the Earth (Ashby et al., 2007). The official science mis- 
sion duration is 1 year with possible extensions to 8 years. 

The accurate measurement of (3 and 7 is one of the sci- 
entific goals also of the proposed mission ASTROD (Astro- 
dynamical Space Test of Relativity using Optical Devices) 
(Ni, 2002; Ni et al., 2004, 2006); the expected accuracy is of 
the order of 10 -7 or better for both (Ni ct al., 2004, 2006). 
The concept of ASTROD is to put two spacecraft in sep- 
arate solar orbits and carry out laser interferometic rang- 
ing with Earth reference stations (e.g. a spacecraft at the 
Earth-Sun L1/L2 points). A simple version of ASTROD, 
ASTROD I, formerly known as Mini-ASTROD, has been 
studied as the first step to ASTROD. ASTROD I employs 
one spacecraft in a solar orbit and carries out interferomet- 
ric ranging and pulse ranging with ground stations. The 
distance of ASTROD I from the Sun varies from 0.5 AU 
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to 1 AU (Ni et al., 2004). The nominal mission lifetime is 
about 1 — 2 year. 

Another very important future mission is LISA (Laser In- 
terferometer Space Antenna) 2 , designed to search for and 
detect gravitational radiation from astronomical sources. 
The LISA concept implies the use of three spacecraft orbit- 
ing the Sun at 1 AU distance in a quasi-equilateral triangle 
formation 20 deg behind the Earth due to a compromise be- 
tween minimising the gravitational disturbances from the 
Earth-Moon system and the communications needs (Pov- 
oleri et al., 2006). The LISA nominal mission is 3.5 years 
with an extended duration of 8.5 years. 

4. The impact of the Kuiper Belt Objects and of 
the minor asteroids 

In order to reach their very ambitious goals, the helio- 
centric orbits of the spacecraft of such missions, and of the 
Earth itself, must be known with high accuracy; any mis- 
modelled or unmodelellcd force acting on them may repre- 
sent a potential systematic bias to be carefully accounted 
for. To date, a dynamical feature which is not yet included 
in the force models of the data reduction softwares used for 
orbit determination purposes is represented by the KBOs. 
Attempts towards the implementation of this goal started 
recently (Pitjeva, 2006; Standish and Williams, 2006). Iorio 
(2007) explicitly worked out the secular effects induced by 
the KBOs on the perihelia of the inner planets pointing out 
that they may affect accurate tests of Einsteinian and post- 
Einstcinian gravity to be performed in that regions of the 
Solar System. Following the approach outlined in (Berto- 
lami and Vieira, 2006; Iorio, 2007), it is possible to analyti- 
cally calculate the cumulative shift in the transverse orbital 
direction T over a time span AP. In the planet's small ec- 
centricity and inclination limit, the Classical KBOs (CK- 
BOs), which follow nearly circular orbits with relatively low 
eccentricities (e < 0.25) and semimajor axes 41 < a < 46 
AU, constituting about 70% of the so far observed popula- 
tion, induce a transverse perturbation 



AT, d AX = -^J^ 7 m * 2+ ^ , (1) 

max "1" -Rmin)-Rmax-^min 

where A is the planetary mean longitude, G is the Newto- 
nian gravitational constant, Mq is the Sun's mass, d and 
e are the planet's semimajor axis and eccentricity, respec- 
tively, tok is the CKBOs' mass, i? max and i? m i n are the 
outer and inner radius, respectively, of the uniform disk 
with which CKBOs have been modelled (Bertolami and 
Vieira, 2006). Eq. (1) comes from the integration over one 
orbital period of the perturbative Gauss equations for the 
variation of the Keplerian elements in terms of which A can 
be expressed (Milani et al., 1987); since for almost circu- 
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lar and equatorial orbits A w M. + SI + u>, where A4 is the 
mean anomaly and Jl is the node, we used 
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where n = ^/ 'GMq/ 'd 3 is the Keplerian mean motion and 
A r is the radial perturbing acceleration (Bertolami and 
Vieira, 2006; Iorio, 2007). 

Eq. (1) and m K = 0.052m e (Iorio. 2007) tell us that 
the nominal values of the transverse shift due to CKBOs 
amounts to -0.10 m for Mercury and -1.01 m for the Earth 
after one year; this corresponds to an average CKBOs accel- 
eration of the order of 10~ 15 m s~ 2 at r « 1 AU. Since the 
uncertainty in the CKBOs' mass is high, Stuk = O.223m,0 
(Iorio, 2007) yields upper bounds of 0.42 m for Mercury 
and 4.33 m for the Earth after one year. In fact, a uni- 
form disk like the one with which the action of the CKBOs 
is modelled induces not only a radial acceleration (Ander- 
son et al., 2002; Nieto, 2005; Bertolami and Vieira, 2006) 
but also a normal one (Owen, 2003). Thus, also the out-of- 
planc portion N of the orbit of a planet/spacecraft is, in 
principle, affected increasing the total orbit bias. However, 
a quantitative evaluation of such a further bias is beyond 
the scope of the present letter. 

Another way to give an order-of-magnitude assessment 
of the bias due to CKBOs on the recovery of the relativistic 
effects of interest, especially useful for BepiColombo since it 
will use, among other things, the orbital motions of Mercury 
and the Earth, is to consider the ratio of the perihelion rates 
induced by CKBOs (Iorio, 2007) to the general relativistic 
ones: it amounts to 6 x 10 -7 for Mercury and 3 x 10 -5 
for the Earth. It may also be interesting to do the same 
for the precessions induced by the solar quadrupole mass 
moment J2, whose measurement at 10 -9 level of accuracy, 
or better, is one of the goals of BepiColombo (Milani et al., 
2002; Ashby et al., 2007): according to Iorio (2007) and by 
assuming J2 = 2 x 10~ 7 , the bias due to the CKBOs is 8 x 
10~ 4 and 1 x 10 _1 for Mercury and the Earth, respectively. 
In Table 1 we release the details. 

The same reasonings can be made for the hidden mass in 
the asteroid belt (Krasinsky et al., 2002) whose action does 
not include the individual perturbations of the 342 major 
asteroids and can be modelled as a massive ring with a con- 
stant mass distribution in the ecliptic plane. By assuming 3 
(Krasinsky et al, 2002) m ring = (5 ± 1) x 10~ 10 M o and 
R k, 2.80 AU, with a 3% uncertainty, we have 



10.7 ± 1 m 
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for a « 1 AU orbit after one year. The ratios of the per- 
ihelion rates induced by the asteroid ring to the general 
relativistic and J2 ones amount to (6.1 ± 1.2) x 10~ 6 and 
(1 ± 0.2) x 10~ 2 for Mercury, and (3 ± 0.6) x 10~ 4 and 
1.25 ± 0.25 for the Earth, as can be noted from Table 1. 



Given the very ambitious goals of the missions examined, 
the figures obtained here may yield some further reasons 
to perform extensive numerical simulations and to enforce 
the process of accurate modelling of the action of KBOs 
and of the minor asteroids in planetary/spacecraft orbit 
determination softwares. In particular, it would be useful 
to estimate the KBOs mass along with the rest of param- 
eters in a global, least-square solution of the Solar System 
dynamics using a massive data set to investigate, among 
other things, the correlation between the parameters of the 
KBO and of the asteroid belt. 
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Table 1 

Nominal values of the secular perihelion precessions lj of Mercury 
and the Earth, in arcseconds per century (" cy — 1 ), induced by general 
relativity (GR), solar quadrupole mass moment (J2), CKBOs and 
the asteroid ring. 

Planet GR J2 CKBOs asteroid ring 
Mercury 42.98047 0.02542 0.00002 0.00026 
Earth 3.83869 0.00084 0.00010 0.00110 



r. 



